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SUMMARY 


The  planfetary  geomagnetic  indices  Kp 


ari^  ap  are  widely  used  in  space 
geophysics,  litey  provide  an  estimate  of  maximum  magnetic  perturbation  within  a 
3-hour  period.  Many  geophysical  properties  are  clearly  related  to  the  indices, 
through  energy  transfer  from  a  common  disturbance  source ,  but  direct  correlation 
is  often  lacking  because  of  poor  matching  between  the  frequency  of  sampling  and 
the  physical  response  functions.  The  index  ap(I)  is  a  simple  accumulation  of 
the  linear  ap  calculated  with  an  attenuation  factor  x  included  to  take 
account  of  natural  temporal  relaxation.  The  case  for  ap(T)  and  the  related 
Kp(t)  is  made  using  applications  to  the  variability  of  the  plasma  environment  in 
the  ionosphere  and  inner  magnetosphere.  These  examples  of  improved  correlation 
suggest  that  tiiue-weighted  integration  jnight  profitably  be  applied  to  other 
indices.  vr  -jLo,  ^ 
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Timc-W'cightcd  Accumulations  ap(x)  and  Kp{x) 


Gordon  L.  Wrenn 


Royal  Aircraft  EMahlishmtnt.  Farnborough,  Hampshire.  England 

The  planetary  geomagnetic  indices  Kp  and  ap  are  widely  used  in  space  geophysics.  They  provide  an 
estimate  of  maximum  magnetic  perturbation  within  a  3-hour  period.  Many  geophysical  properties  arc 
clearly  related  to  the  indices,  through  energy  transfer  from  a  common  disturbance  source,  but  direct 
correlation  is  often  lacking  iKcause  of  poor  matching  between  the  frequency  of  sampling  and  the  physical 
response  functions.  The  index  up(t(  is  a  simple  accumulation  of  the  linear  ap  calculated  with  an  attenu¬ 
ation  factor  r  included  to  lake  account  of  natural  temporal  relaxation.  The  case  for  apfr)  and  the  related 
Kp(T]  is  made  using  applications  to  the  variability  of  the  plasma  environment  in  the  ionosphere  and 
inner  magnetosphere.  These  examples  of  improved  correlation  suggest  that  time-weighted  integration 
might  profitably  be  applied  to  other  indices. 


1.  Introduction 

Geomagnetic  indices  are  indispensable  to  much  of  space 
geophysics,  but  they  do  have  severe  inherent  limitations  which 
are  sometimes  ignored.  In  advancing  a  new  treatment,  it  is 
important  to  appreciate  the  underlying  constraints  which  will 
be  carried  through.  Rnsinker  [1972]  concludes  that  “indices 
are  of  errai  use  in  statistical  studies  of  solar-terrestrial  interac¬ 
tions.  where  qualitative  correlations  between  different  parame¬ 
ters  are  being  sought"  but  "it  is  not  advisable  to  use  indices  in 
a  study  of  individual  events,  particularly  where  one  wishes  to 
rule  out  the  possibility  of  substorm  activity.”  He  also  notes 
that  "in  general,  indices  are  only  capable  of  defining  the  lower 
limit  of  geomagnetic  activity  at  any  given  time." 

The  3-hour  ap  index  is  a  linear  scaling  of  Kp  which  is 
derived  from  a  statistical  composite  of  the  variations  at  a 
selected  group  of  subauroral  zone  stations  IBartels  and  Veld- 
kump.  1949],  The  basic  input  data  are  the  differences  between 
absolute  maximum  and  absolute  minimum  for  the  three  com¬ 
ponents  of  the  magnetic  flux  density,  measured  in  nanoteslas, 
in  each  3-hour  interval  [Riwioker,  1972],  The  computation, 
involving  local  K  and  Ks  values  and  observatory  codings  ap¬ 
propriate  to  season  and  local  lime,  is  complex  [Mayaud, 
1980]  but  the  consequence  is  that  ap  can  be  numerically  relat¬ 
ed  to  the  magnitude  of  the  disturbance  at  a  "standard"  mid- 
latitude  station.  The  scale  of  ap  is  0  to  400  (that  of  Kp  is  Oo  to 
9o|.  and  each  unit  is  approximately  equivalent  to  a  flux  den¬ 
sity  change  of  2  nT.  There  are  numerous  problems  with  de¬ 
tailed  scaling  of  magnetograms,  but  a  potential  limitation  of 
ap  and  Kp  is  introduced  by  the  restricted  and  nonuniform 
distribution  of  observatories  around  the  Earth  with  respect  to 
longitude,  local  time,  and  geomagnetic  latitude  within  the 
chosen  band  between  47  and  63  [Michel,  1964], 

The  merit  of  the  ap  index  lies  in  its  widespread  availability 
and  the  continuity  which  now  spans  more  than  50  years. 
There  now  exist  a  great  many  empirical  relationships  between 
ap  (Kp)  and  other  geophysical  parameters;  this  considerably 
extends  the  uttefulness  of  the  index.  Other  indices  (e.g.,  Dst, 
AV.  A  I,.  AE.  r.)  and  direct  particle  and  field  measurements  are 
available  and  may  be  more  appropriate  for  particular  studies. 
A  number  of  new  indices  have  been  proposed,  but  it  .s  a 
matter  of  some  controversy  whether  they  are  needed  or 
whether  they  could  supersede  established  ones  [Mayaud, 
1980], 

Published  in  1987  by  the  American  Geophysical  Union. 

Paper  number  7A906I. 


The  index  opfr)  is  a  time  series  accumulation  of  ap  with 
exponential  smoothing.  It  is  an  attempt  to  produce  an  index 
which  reflects  an  integration  of  geomagnetic  activity  over  a 
number  of  3-hour  intervals,  giving  more  weight  to  the  recent 
past  and  relatively  less  to  measurements  from  earlier  times.  It 
has  been  used,  with  some  success,  to  overcome  certain  inad¬ 
equacies  in  the  common  indices  by  the  introduction  of  the 
integrating  and  smoothing  properties  coupled  with  a  matched 
persistence  factor.  The  classification  of  data  as  "quiet"  or  “dis¬ 
turbed"  must  be  arbitrary,  but  it  has  been  variously  based 
upon  index  values,  sums,  or  averages  for  selected  periods  prior 
to  the  events  of  interest.  Wrenn  et  al.  [1987]  introduced  ap<r) 
to  control  the  binning  of  values,  first  to  establish  valid 
quiet  time  references  and  then  to  extract  main  phase  storm 
patterns  from  their  data.  Here  it  is  described  in  detail  with 
additional  examples  to  demonstrate  that  it  can  also  be 
matched  to  studies  of  the  inner  magnetosphere  and  related 
plasma  effects. 

2.  Derivation  of  3-Hourly  a,  ,r)  and  Kp(t) 

Values  of  ap  have  been  tabulated  for  all  years  since  1932 
and  are  now  available  on  computer  input  media.  Any  integra¬ 
tion  is  limited  by  the  3-hour  resolution  and  cannot,  therefore, 
be  usefully  applied  to  events  with  time  scales  of  le.ss  than 
several  hours.  A  3-hour  attenuation  multiplier  r  is  introduced 
(0  <  t  <  1),  and  aplr)  defined  by  a  geometric  progression  to  be 

ap{r)  =  (I  -  i)[op  -E  (T)ap. ,  +  (y‘)ap .  ^  -E  ■] 

where  ap  _ , .  ap  _  j  are  ap  values  for  —  3  hours,  —  6  hours,  etc. 

The  factor  ( I  —  r)  normalizes  the  summation  because 

I  -E  r  E  -E  T^ -E  ■  =(1  -r)-' 

The  repeated  summing  of  infinite  series  can  be  avoided  by 
setting  up  a  recursion  loop  and  simply  stepping  sequentially 
through  the  table  of  ap  values,  applying  the  multipler  as 

I - -  op(r)  =  [I  -  r]  np  -E  [T]ap  ,(t) - s 


An  initial  value  of  ap  _  ,(t)  is  required.  If  this  is  unavailable, 
then  a  dummy  number  (say,  6)  may  he  input,  in  which  case  the 
output  for  the  first  few  days  should  be  disregarded. 

For  every  value  of  r.  ap(r)  falls  within  the  same  range  of 
values  as  ap,  0  to  4(X>.  Note  that  ap(0)  equals  the  single  un¬ 
smoothed  3-hnur  value  of  ap. 
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Fig.  I.  V«ri<lions  ofy^Fj  and  ap(x)  for  November  17-30.  1971. 
The  lop  panel  shows  In  N/Af,  for  Argentine  Islands.  Electron  con¬ 
centration  S  is  computed  from  hourly  values,  and  Onset)  is 
interpolated  from  the  monthly  medians,  for  November  (curve)  and 
December  (points):  units  are  cm*’.  Other  panels  show  ap(x}  for 
t  -  0.0. 0.S,  0.75,  and  0.9. 

In  order  tq  demonstrate  the  form  of  ap{x).  Figure  I  plots 
four  profiles  for  the  second  half  of  November  1971  using 
values  of  T  0.0, 0.S,  0.7S,  and  0.9.  The  effect  of  increasing  t  is 
clear,  as  the  short-term  variations  are  smoothed  out  in  the 
integrating  operation.  The  curves  are  inverted  to  ease  the 
comparison  with  the  changes  in  F  region  electron  density  pre¬ 
sented  in  the  same  figure  and  discussed  in  the  next  section. 

The  choice  of  multiplier  r  is  fundamental  to  the  procedure 
and  depends  upon  the  characteristic  response-recovery  time  of 
the  effect  being  investigated.  This  is  best  illustrated  by  taking  a 
delta  function  in  ap  and  comparing  the  subsequent  decay  of 
ap(x)  for  different  values  of  t.  Figure  2  presents  such  a  plot  for 
a  10-day  period  using  t  0.99, 0.975, 0.95, 0.90, 0.75,  and  0.50. 
The  times  required  for  a  l/<  decay  are  12.4  days,  4.9  days,  2.4 
days,  1.2  days,  10.4  hours,  and  4.3  hours,  respectively,  ft  is 
useful  to  define  persistence  as  3A1  —  r)  hours  as  an  approxi¬ 
mate  measure  of  these  timet  (see  Table  I). 

ACpfr)  can  be  determined  from  ap(x)  using  the  standard  con¬ 
version  table  with  appropriate  bisectors.  Kp{x)  does  permit  a 
direct  comparison  of  accumulated  values  with  those  used  pre¬ 
viously  in  many  applications.  Ap  is  a  daily  average  of  ap, 
being  a  simple  mean  of  the  eight  3*hourly  ap  values.  J.  H. 
Allen  (personal  communication,  1985)  recognized  a  problem  of 
using  Ap  to  characterise  magnetic  storms  which  start  well  into 
the  UT  day  and  introduced  Ap*  as  the  largest  value  computed 
from  successive  eight-point  running  means  of  ap.  The  author 


DAYS 


Fig.  2.  Decay  of  ap(x)  for  x  -  0.5. 0.75. 0.9. 0.95.  0.975.  and  0.99.  The 
contribution  of  any  single  value  of  ap  reduces  with  time. 

knows  of  no  other  attempt  to  systematically  generate  a  retro- 
gressively  weighted  sampling  of  geophysical  indices,  but  there 
have,  of  course,  been  numerous  searches  for  lagged  relation¬ 
ships. 

Although  the  scale  of  ap{x)  is  the  same  as  that  o'  ap,  the 
distribution  of  values  is  changed  by  the  smoothing,  extremes 
being  eliminated.  This  is  shown  by  an  examination  of  ex¬ 
tremely  disturbed  times.  The  most  active  period  in  recent 
years  occurred  in  February  1986,  when  ap  registered  as  fol¬ 
lows: 


Feb.  7 

Feb.  8 

Feb.  9 

15 

132 

300 

22 

132 

179 

18 

154 

56 

67 

111 

39 

67 

207 

48 

179 

179 

94 

III 

400 

56 

179 

300 

27 

Determination  of  ap(r)  leads  to  the  following  maxima: 

op(0.50)  =  297 

[Xp(0.50)  = 

9-] 

flp(0.75)  =  252 

[Kp(0.75)  = 

8  +  ] 

ap(0.90)  =  167 

[/Cp(0.90)  - 

8-] 

ap^0.95)  =  108 

[Kp(0.95)  - 

7-] 

These  figures  can  be  compared  with  those  for  the  previous 
superactive  period  in  August  1972,  when  the  maxima  were 

ap(0.50)  -  240 

[Kp(0.50)  - 

8  +  ] 

op(0.75)  =  204 

[Ap(0  75)  - 

8o] 

ap(0.90)  -  143 

[/Cp(0.90)  - 

7  +  ] 

ap(0.95)  =  1 00 

[Ap(0.95)  = 

6  +  ] 

TABLE  I.  Frequency  Distribution  of  ap(r),  1974-1984 


T 

Persistence, 

hours 

apix) 

Peak 

10% 

25% 

75% 

OO”/!. 

0 

3 

9 

2.8 

4.9 

17.5 

33.3 

0.5 

6 

5.5 

3.8 

5.8 

19.7 

32.8 

0.75 

12 

6.0 

4.4 

6.7 

19.9 

31.8 

0.9 

30 

7.5 

5.6 

8.0 

20.5 

29.7 

0.95 

60 

10.0 

6.8 

9.2 

19.7 

27.8 

0.975 

120 

11.5 

8.1 

10.4 

19.4 

25.9 

0.99 

300 

13.5 

9.9 

11.4 

19.0 

23.4 
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Bartels  [1963]  discussed  the  variations  of  Kp  and  up  over 
three  solar  cycles  and  defined  frequency  distributions  for  se¬ 
lected  epochs.  With  apix),  the  frequency  distribution  is  a  func¬ 
tion  of  r,  and  this  has  been  examined  for  one  unexceptional 
solar  cycle,  1974  to  1984  (32.144  intervals).  The  following  divi¬ 
sions  were  set:  very  quiet,  ;  quiet,  10-25%;  disturbed, 

75-90%;  and  very  disturbed.  90  100%.  Half  the  distribution, 
25-75%,  could  then  be  classed  as  normal  or  moderately  dis¬ 
turbed.  Employing  a  step  width  of  0.5,  the  frequencies  were 
plotted  for  the  selected  values  of  i.  Table  1  lists  the  interpo¬ 
lated  values  of  ap(z),  corresponding  to  the  frequency  peak  and 
the  prescribed  percentile  boundaries,  for  the  same  set  of  r 
values  as  was  used  in  Figure  2. 

The  75%  disturbed  threshold  remains  remarkably  stable 
(Kp  =  3-1-)  but  the  frequency  distribution  does  change  with 
increasing  persistence,  and  this  must  be  taken  into  account 
when  comparing  results  from  analyses  which  use  diflereni 
values  of  T. 

Monthly  listings  of  Kp  are  distributed  routinely  by  the  In- 
stitut  fiir  Geophysik.  Gottingen,  Federal  Republic  of  Ger¬ 
many,  and  are  readily  available  from  World  Data  Centres. 
The  derivation  of  ap(r)  is  very  simply  accomplished  in  a  com¬ 
puter  code  or  within  an  instruction  sequence  on  a  program¬ 
mable  calculator. 

3.  f  RECiioN  Storms 

The  spur  for  the  development  of  ap(T)  came  from  a  study  of 
ionospheric  storms  and  depictions  in  the  electron  con¬ 
centration  at  the  peak  of  the  F  layer,  i.e..  negative  storm  ef¬ 
fects.  From  ionosonde  data,  it  is  found  that  the  critical  fre¬ 
quency  /,Fj,  measured  at  mid-latitude  stations,  can  decrease 
considerably  with  respect  to  undisturbed  values  for  several 
days  at  times  of  geomagnetic  activity  [Rishbeth.  1975].  This  is 
illustrated  by  the  top  panel  of  Figure  1.  where  hourly  values  of 
/„Fj  are  converted  to  electron  concentrations  N  and  plotted  as 
In  (/V/A/q),  where  is  an  undisturbed  value  derived  from  the 
appropriate  monthly  medians.  The  data  from  Argentine  Is¬ 
lands  (54°S  geotfiagnetic)  refer  to  November  1971.  The  diurnal 
median  profiles  of  for  November  and  December  are  in¬ 
cluded  in  the  inset.  The  depletions,  strongest  during  the  morn¬ 
ing  hours,  are  effectively  confined  to  the  period  of  the  20th  to 
26th  when  the  high  geomagnetic  disturbance  is  reflected  by  up. 
This  qualitative  correlation  has  been  demonstrated  in  numer¬ 
ous  studies  of  ionospheric  storms.  Treading  carefully  through 
many  similar  “events.”  it  is  clear  that  high  values  of  ap  do 
coincide  with  Fj  storm  effects,  but  the  latter  do  not  persist 
unless  the  disturbance  continues,  i.e.,  unless  ap  remains  high. 
The  extent  of  the  storm  effect  also  depends  more  on  the 
average  value  of  ap  through  the  event  rather  than  the  maxi¬ 
mum  which  might  be  limited  to  a  single  3-hour  interval.  Sig¬ 
nificant  storm  effects  are  observed  when  ap  stays  at  a  moder¬ 
ate  value  for  several  days  without  a  sudden  commencement  or 
a  high  peak.  The  diurnal  variations  were  clear  and  their  en¬ 
velope  depended  upon  latitutde  and  season  as  well  as  activity, 
but  these  effects  could  not  be  investigated  without  a  suitable 
“activity  index”;  ap  and  Dsi  proved  unsatisfactory. 

The  logic  of  trying  to  integrate  ap  is  obvious  even  though 
any  single  value  reflects  the  peak  deviation  in  a  3-hour  period. 
The  physical  processes  producing  the  storm  effects  are  not 
understood,  but  some  modification  of  the  circulation  of  global 
thermospheric  winds  is  thought  to  be  responsible  [Volland, 
1983;  Sojka  and  Schunk,  1984].  This  relates  directly  to  the 
atmospheric  energy  input  which  could  be  associated  with  a 


suitable  integration  or  accumulation  of  ap.  Any  effect  of  an 
energy  input  will  last  for  a  limited  time,  and  this  natural  relax¬ 
ation  can  be  built  into  the  accumulation. 

Plots  using  the  formats  of  Figure  1  were  produced  for  five 
stations  for  8  years.  The  negative  storm  effects  are  complex, 
but  consistent  patterns  have  been  identified  by  sifting  the  data 
with  respect  to  ap(z).  Again  it  must  be  stressed  that  a  one  to 
one  correlation  on  single  events  is  beyond  the  scope  of  the 
indices.  Initially,  t  =  0.9  was  chosen  because  this  well  matches 
the  common  recovery  response  of  Dst,  but  it  became  clear 
from  the  whole  survey  that  the  F  region  generally  recovers 
faster  than  this  and  the  appropriate  value  of  t  is  closer  to  0.75. 
The  curves  displayed  in  Figure  I  do  support  this  conclusion, 
but  it  may  be  dilficult  to  see  without  producing  overlays. 

The  power  of  the  new  index  is  demonstrated  by  treating  the 
8  years  of  Argentine  Islands  data,  selecting  summer  months 
(November-January),  and  binning  all  the  computed  In  (N/N„) 
values  in  terms  of  hour  and  ap(0.TS)  ranges  0-7,  7-20,  20-32, 
and  32-400.  Figure  3  presents  the  results  and  establishes  clear 
diurnal  depletion  patterns  for  disturbed  (20-32)  and  very  dis¬ 
turbed  (>32)  conditions.  It  is  interesting  that  the  monthly 
median  reference  values  are  associated  with  “moderate”  dis¬ 
turbance  (7-20)  so  that  quiet  and  very  quiet  intervals  (0-7) 
generate  apparent  positive  effects.  This  is  a  neat  demon¬ 
stration  of  the  inadequacy  of  monthly  medians  as  a  quiet  time 
reference. 

4.  Plasma  Environment  at  6.6  R^: 

Spacecraft  Anomalies 

Plasma  boundaries  observed  at  geosynchronous  orbit  (6.6 
r,  )  arc  clearly  subject  to  geomagnetic  control.  Cold  plasma, 
populating  the  plasmasphere,  is  convected  away  in  disturbed 
periods;  the  refilling  then  takes  several  days.  Protons  and  elec¬ 
trons,  with  characteristic  energies  of  less  than  I  eV,  are  gener¬ 
ally  observed  in  the  postnoon-evening  sector  but  are  fre¬ 
quently  found  at  other  local  times.  Refilling  of  the  flux  tubes 
takes  many  days,  but  Wrenn  el  ai  [1984]  have  supported  the 
notion  of  filling  from  the  center;  this  would  suggest  an  equa¬ 
torial  persistence  more  like  30  hours. 

GEOS  2  measured  cold  plasma  concentration  at  geostation¬ 
ary  orbit  between  August  1978  and  July  1980.  The  marked 


STATION  TIME  (hr) 


Fig.  3.  Diurnal  patterns  of  main  phaK  storm  effects.  Hourly 
mean  values  of  In  (JV/Af,)  from  Argentine  Islands  (summer,  1971- 
1978)  are  plotted  within  selected  ranges  of  ap(0.75).  The  average  num¬ 
bers  of  samples  per  point  were  as  follows;  ap(0.7S)-  32-400,  50; 
ap(0.75)  =  20-32, 103;  opfO.TJ)  -  7-20, 349;  apiO.JS)  -  0-7,  203. 
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day  to  day  variability  is  clearly  related  to  geomagnetic  ac¬ 
tivity,  but  neither  ap  nor  Ap  has  proved  suitable  for  empirical 
modeling.  However,  averaging  the  data  with  respect  to  ap(t), 
does  give  meaningful  results.  Figure  4  uses  t  =  0.9  to  split  the 
data  into  (from  top  to  bottom)  very  disturbed,  disturbed, 
normal,  quiet,  and  very  quiet  and  plots  average  concentration 
(in  cm  ~  profiles  as  a  function  of  local  time.  The  proton  data 
are  maximum  values  detected  within  hourly  intervals.  The 
curves  neatly  describe  the  broad  features  of  the  equatorial 
plasmasphere:  it  extends  beyond  6.6  at  all  local  times  after 
a  period  of  sustained  quiet  [Sojka  and  Wrenn,  1985]  and  nor¬ 
mally  exhibits  a  bulge  near  1800  LT,  but  concentrations  peak 
at  earlier  local  times  during  active  spells.  It  has  been  possible 
to  demonstrate  the  worth  of  using  t  =  0.9  rather  than  t  =  0 
for  fitting  the  data,  but  the  set  is  too  small  to  yield  a  statis¬ 
tically  significant  “best  value'. 

Geosynchronous  spacecraft  frequently  suffer  from  surface 
charging  effects;  the  abundance  of  cold  plasma  is  important 
because  it  prevents  hazardous  charging,  but  in  the  past,  it  has 
been  consistently  omitted  from  environment  characterisations 
IGarrett  el  al.,  1981].  Hot  plasma  sheet  electrons  are  usually 
confined  to  the  nightside,  but  fluxes  increase  in  both  intensity 
and  duration  as  geomagnetic  activity  rises;  spacecraft 
charging  effects  and  operational  anomalies  then  become  more 


LOCAL  TIME 

Fig  4.  Cold  plasma  concentration  at  6i6  Rf.  Average  daily  pro¬ 
files.  determined  from  GEOS-2  measurements  between  August  1978 
and  July  19(0  are  shown  for  five  ranges  of  apiO.9)  from  very  quiet 
I  <  S.5)  to  very  disturbed  ( >  30). 
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Fig.  5.  Mcieosat  2  anomaly  analysis  using  Ap  and  ap(0.9S). 
Forty-seven  radiometer  scan  anomalies  occurred  between  September 
1981  and  September  I98S.  Their  distribution  with  respect  to  geomag¬ 
netic  activity  is  shown  (a)  as  a  percentage  of  days  in  each  Ap  range 
(1491  days)  and  (b)  as  a  percentage  of  intervals  in  each  3-hourly 
np(0.95)  range  (11,928  intervals).  The  vertical  scales  are  equivalent 
when  the  sampling  factor  of  8  is  taken  into  account. 

frequent.  This  phenomenon  provides  a  significance  test  for 
ap(t)  with  a  value  of  r  close  to  1. 

Between  September  1981  and  September  1985,  Meteosat  2 
suffered  47  radiometer  scan  anomalies  which  interrupted  the 
continuous  Earth  imaging  function.  Were  these  anomalies  re¬ 
lated  to  geomagnetic  conditions?  The  rather  uniform  distri¬ 
bution  in  local  time  precludes  any  direct  triggering  by  plasma 
sheet  electron  fluxes,  but  more  indirect  mechanisms  cannot  be 
ruled  out.  Figure  S  shows  the  result  of  ar  analysis  undertaken 
with  respect  to  the  daily  Ap  index.  The  1491  days  were  binned 
with  width  5,  and  for  each  bin,  the  fraction  of  days  exhibiting 
an  anomaly  was  calculated.  Figure  Sa  suggests  a  bias  towards 
the  larger  Ap  values  but  it  is  not  very  convincing.  Johnstone  et 
at.  [1985]  concluded  that  the  “anomalies  tend  to  occur  after 
an  extended  period  ( >  4  days)  of  more  than  average  magnetic 
activity.'  Repeating  the  analysis  using  ap(0.95)  for  the  1 1,928 
3-hour  intervals  gives  Figure  5b,  which  provides  clear  evidence 
that  the  number  of  anomalies  does  depend  upon  the  level  of 
activity  appertaining  to  the  previous  60  hours,  although  there 
is  no  serious  attempt  to  determine  the  persistence  here. 

5.  Discussion 

The  objective  of  this  brief  paper  is  simply  to  introduce  the 
concept  of  ap(T)  and  to  demonstrate  its  potential  as  an  aid  in 
working  with  the  common  planetary  indices.  A  detailed  analy¬ 
sis  of  its  application  with  F  region  storms  is  given  by  Wrenn  ei 
al.  [1987].  The  examples  from  geosynchronous  orbit  were  se¬ 
lected  because  the  response  times  appeared  to  be  much  longer 
than  1  day.  In  such  cases,  Ap  or  Y,  could  be  adequate  if 
the  tie  to  the  UT  day  is  removed  h  upjr)  permits  an  opti¬ 
mization  with  the  selection  of  r,  and  it  allows  for  the  retention 
of  diurnal  variations  in  the  analyses. 

When  applying  a  lime  dependent  filter  to  data,  it  is  essential 
to  consider  whether  the  product  is  distorted;  characteristic 
frequencies  of  the  filter  should  not  intrude.  Here  the  ap  input 
data  are  square  pulses,  sequential  averaging  merely  serves  to 
attenutate  the  basic  3-hour  filter.  The  normalized  series  is  a 
discrete  first-order  one  of  a  simple  Markoff  type;  it  has  been 
used  widely  in  autocorrelation  and  auloprojective  analysis 
IKendall,  1973,  p.  1 18]. 
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The  Fourier  transform  of  this  RC  filter  can  be  defined  as 


ill- 
2>t/(l  -  r 


where  T  =  3  hours  in  this  case. 

As  T  increases  from  0  to  1  the  transform  changes  smoothly 
from  T  sine  (n/T)  to  X/inf.  Given  the  nature  of  ap,  frequencies 
greater  than  1/T  are  only  of  academic  interest.  Within  the 
range  of  application,  say  O.OI/T  to  1/T,  the  shape  of  the  spec¬ 
tral  window  is  changed  as  t  increases,  but  no  “side  lobes"  are 
introduced.  The  important  natural  triggering  frequency  here  is 
0.125/T  (i.e.,  24-hour  period). 

Indices  are  an  essential  ingredient  of  empirical  modeling, 
but  they  often  fall  short  of  being  an  effective  aid  to  the  under¬ 
standing  of  relevant  physics.  If  the  frequency  of  the  index  is 
small  compared  with  those  characterizing  the  phenomenon 
being  studied  then  such  an  index  is  clearly  inappropriate.  Ac¬ 
cordingly,  the  need  for  I -min  AE  has  been  appreciated  in  the 
study  of  substorms.  When  the  frequency  of  the  index  is  com¬ 
parable  or  relatively  high,  then  it  is  still  important  to  achieve  a 
good  match  between  the  sampling  of  the  index  and  the  time 
constants  involved  in  the  process.  The  type  of  accumulation 
described  here  for  ap  might  well  be  applied  to  AE.  AV,  or  AL 
for  the  investigation  of  subslorm  responses  in  the  1-  to  100- 
min  range. 

Numerous  parameters  have  been  correlated  with  ap  or  Kp 
using  a  lag  time,  but  it  would  now  be  instructive  to  include  r 
as  a  variable  in  the  regression  analyses  in  attempts  to  establish 
appropriate  response  times:  these  could  provide  valuable 
pointers  to  the  identification  of  the  controlling  physical  pro¬ 
cesses.  It  is  possible  that  the  coarseness  or  nonlinearity  in  ap 
might  preclude  a  meaningful  result,  but  given  a  very  large  data 
base,  the  broad-brush  statistical  treatment  is  surprisingly 
powerful. 

Relationships  between  different  indices  are  of  special  in¬ 
terest.  Dsi  monitors  the  storm  time  enhancement  and  subse¬ 
quent  decay  of  the  ring  current;  the  recovery  phase  typically 
lasts  for  3  days.  An  attempt  to  correlate  Ds(  and  apfr)  could 
throw  more  light  on  the  real  significance  of  both  indices  and 
the  relevant  processes  of  particle  injection,  trapping,  and  pre¬ 
cipitation. 

It  is  hoped  that  the  few  examples  presented  here  will  con¬ 
vince  space  scientists  that  there  can  be  significant  merit  in 
employing  a  simple  accumulation  procedure  to  extend  the 
usefulness  of  the  readily  available  geomagnetic  indices:  ap  is  a 
crude  index  of  geomagnetic  activity  and  one  must  be  aware  of 
its  limitations:  the  3-hour  resolution  is  clearly  inadequate  for 
many  studies,  and  better  indices  may  be  available.  Even  so. 
the  introduction  of  weighted  running  averages  does  offer  some 
advantage  in  the  search  for  correlations  with  time-varying  fea¬ 


tures  of  the  ionosphere  thermosphere,  the  plasmasphere,  and 
the  inner  magnetosphere,  and  upfrl  represents  a  useful  tool  in 
the  building  of  a  better  understanding  of  the  physics  of  some 
complex  geospace  systems  from  raw  materials  supplied  by 
routine  monitoring. 

Aiknowleiigmenis  This  work  was  started  at  the  University  of 
Southampton  with  the  support  of  NERC  fprojecl  GST.OZI  12).  Henry 
Rishbelh  and  Alan  Rodger  are  thanked  for  inspiring  the  research, 
encouraging  the  ideas  which  are  formulated  here,  and  suggesting  re¬ 
finements  to  the  script.  The  MuUard  Space  Science  Laboratory  pro¬ 
cessed  the  GEOS  and  Meteosat  data:  Alan  Johnstone's  work  with  the 
latter  is  gratefully  acknowledged.  Mike  Hapgood  of  the  World  Data 
Centre  Cl  supplied  the  ionosonde  scalings  and  geomagnetic  indices 
on  tape  and  floppy  disc. 

The  Editor  thanks  H.  W.  Kroehl  and  M.  Sugiura  for  their  assist¬ 
ance  in  evaluating  this  paper. 

References 

Bartels.  J..  Discussion  of  time-variations  of  geomagnetic  activity,  indi¬ 
ces  Kp  and  ap,  1932-1961.  Ann.  Gtopkys.,  19.  1-20.  1963. 

Bartels,  J..  and  J.  Veidkamp.  Geomagnetic  and  solar  data,  J.  Geophvs. 
Res..  54.  295-299.  1949. 

Garrett.  H.  B..  D.  C.  Schwanfc.  and  S.  E.  DeForest,  A  statistical 
analysis  of  the  low-energy  geosynchronous  plasma  environment.  II. 
Ions!  PXanel.  Space  Sci..  29.  1045-1060,  1981, 

Johnstone.  A.  D.,  G.  L.  Wrenn.  H.  E.  Huckle.  and  R.  F.  Scott,  Me- 
leosat  F2  spacecraft  charging  monitors.  Final  report,  Eur.  Space 
Agency  contracts  47I5'8I;F  CG  and  5911/84/F  CG,  269,  pp.,  Mul- 
lard  Space  Sci.  Lab..  Univ.  Coll.  London.  1985 
Kendall.  M.  G..  Time-Series.  197  pp..  Griffin.  London.  1973. 

Mayaud.  P.  N..  Derivation.  Meaning,  anil  Use  oj  Geomagnetic  Indices, 
Geophvs.  Monogr.  Ser..  vo).  22.  154  pp..  AGU.  Washington,  D.  C.. 
1980.  ' 

Michel.  F  C..  Kp  as  a  planetary  index.  J.  Geophvs.  Res.,  69.  4182- 
4183.  1964. 

Rishbelh.  H..  F-region  storms  and  thermospheric  circulation,  J. 

Atmos.  Terr.  Rhys..  S7.  1055-1064.  1975. 

Rostoker.  G..  Geomagnetic  indices.  Rev.  Geophys..  10. 935-950.  1972. 
Sojka.  J.  J..  and  R.  W.  Schunk.  A  theoretical  F  region  study  of  ion 
composition  and  temperature  variations  in  response  to  mag- 
neiospheric  storm  inputs.  J  Geophys.  Res..  89,  2348-2358.  1984. 
Sojka.  J.  J..  and  G.  L.  Wrenn.  Refilling  of  geosynchronous  flux  tubes 
as  observed  at  the  equator  bv  GEOS  2.  J.  Geophys.  Res.,  90,  6379- 
6385.  1985. 

Volland.  H..  Dynamics  of  the  disturbed  ionosphere.  Space  Sci.  Rev., 
34.  327-335.  1983 

Wrenn.  G  L..  J.  J.  Sojka.  and  J.  F  E.  Johnson.  Thermal  protons  in 
the  morning  magnetosphere:  Filling  and  healing  near  the  equa¬ 
torial  plasmapause.  Planet.  Space  Sci..  32.  351-363.  1984. 

Wrenn.  G.  L..  A.  S.  Rodger,  and  H.  Rishbelh,  Geomagnetic  storms  in 
the  Antarctic  F-region.  I.  Diurnal  and  seasonal  patterns  for  main 
phase  effects.  J.  Atmos.  Terr  Phys..  49.  901-913.  1987. 


G  L.  Wrenn.  Space  Department.  Royal  Aircraft  Establishment. 
Farnborough.  Hampshire.  GUI4  6TD.  England. 

•  Received  April  10.  1987: 
revised  June  5.  1987; 
accepted  June  9. 1987.) 


I  The  planetary  geomagnetic  indices  Kp  and  ap  are  widely  used  in  space  j 

I  geophysics.  They  provide  an  estimate  ot  maximum  magnetic  perturbation  within  a 
3-hour  period.  Many  geophysical  properties  arc  clearly  related  to  the  indices, 
through  energy  transfer  from  a  common  disturbance  source,  but  direct  correlation  is 
often  lacking  because  of  poor  matching  between  the  frequency  of  sampling  and  the 
physical  response  functions.  The  index  ap(T)  is  a  simple  accumulation  of  the 
linear  ap  calculated  with  an  attenuation  factor  t  included  to  take  account  of 
natural  temporal  relaxation.  The  case  for  ap(T)  and  the  related  Kp(T)  is  made 
using  applications  to  the  variability  of  the  plasma  environment  In  the  Ionosphere 
and  inner  magnetosphere.  These  examples  of  improved  correlation  suggest  that 
time-weighted  integration  might  profitably  be  applied  to  other  indices. 


